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Abstract — Time synchronization is important for a variety of 
applications in wireless sensor networks including scheduling 
communication resources, coordinating sensor wake/sleep cycles, 
and aligning signals for distributed transmission/reception. This 
paper describes a non-hierarchical approach to time synchroniza- 
tion in wireless sensor networks that has low overhead and can be 
implemented at the physical and/or MAC layers. Unlike most of 
the prior approaches, the approach described in this paper allows 
all nodes to use exactly the same distributed algorithm and does 
not require local averaging of measurements from other nodes. 
Analytical results show that the non-hierarchical approach can 
provide monotonic expected convergence of both drifts and offsets 
under broad conditions on the network topology and local clock 
update stepsize. Numerical results are also presented verifying 
the analysis under two particular network topologies. 

Index Terms — synchronization, timing, wireless sensor net- 
works, consensus clock 

I. Introduction 

Synchronization is the process of establishing a common 
notion of time among two or more entities. In the context of 
wired and wireless communication networks, synchronization 
enables coordination among the nodes in the network and can 
facilitate scheduling of communication resources, interference 
avoidance, event detection/ordering, data fusion, and coordi- 
nated wake/sleep cycles [1|. Precise synchronization, to the 
order of a fraction of a carrier period, can also enable effi- 
cient distributed transmission schemes such as retrodirective 
distributed beamforming 0. 

Standardized protocols for synchronizing devices in a net- 
work include Network Time Protocol (NTP) 1 3 1 and Precision 
Time Protocol (PTP, also known as IEEE 1588) [4J. Both NTP 
and PTP are hierarchical synchronization systems: devices in 
the network are assigned to a class or stratum and those 
with the lowest stratum number are assumed to be perfectly 
synchronized with Coordinated Universal Time (UTC). Nodes 
with higher stratum numbers synchronize their clocks via 
TCP/IP messages with nodes having lower stratum numbers. 
While NTP is currently responsible for synchronizing the 
clocks of most of the devices connected to the internet, it (and 
PTP) are generally considered too cumbersome for networks 
of inexpensive computationally-constrained sensor nodes. 

Several synchronization protocols have been developed to 
address some of the shortcomings of NTP in a sensor- 
network scenario. Protocols based on bidirectional messaging 
between nodes in adjacent levels of a hierarchical tree include 
the timing-sync protocol for sensor networks (TSPN) 0, 
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lightweight tree-based synchronization (LTS) |6|, and tiny- 
sync/mini-sync [7 |. While the details of these protocols differ, 
they are all based on the establishment of a hierarchical tree 
structure and the goal is to synchronize every node in the 
network to the root node. The use of bidirectional messages 
allows for disambiguation of propagation delay from clock 
offset, making these protocols potentially more accurate than 
synchronization schemes based on unidirectional messages. 

The flooding time synchronization protocol (FTSP) [HO, 
reference broadcast synchronization (RBS) (9), and Cesium- 
Spray [10 1 are among the most widely cited examples of 
synchronization protocols for sensor networks based on unidi- 
rectional messaging. While the use of unidirectional messages 
potentially reduces overhead, it prevents propagation delay 
compensation and limits the attainable synchronization accu- 
racy. FTSP specifies a hierarchical structure with timestamped 
messages flowing from the root node to nodes in lower levels 
of the tree. RBS and CesiumSpray are based on broadcast 
beacon transmissions followed by the exchange of timestamps 
among receivers in the broadcast range of the beacon node. 
In networks with more than one beacon node, RBS gateway 
nodes are specified to transform timestamps between broadcast 
domains. While RBS and Cesium Spray are non-hierarchical, 
they still require some network structure (and corresponding 
overhead) in the establishment of special beacon nodes, the 
establishment of the sets of client nodes within the broadcast 
range of each beacon node, and the exchange of timestamps 
among the set of client nodes as pointed out in ifTTl . 

A handful of synchronization techniques based on consen- 
sus and diffusion have also recently been proposed. While 
some of these techniques require network hierarchy, e.g. \12\, 
iTDl . others are non-hierarchical and use local averaging 
of clock values reported by neighboring nodes to achieve 
consensus (TJ, lfT3l - lfl31 . All of these techniques are based 
on MAC layer timestamping. Physical layer synchronization 
techniques have also recently been reported including carrier 
synchronization for coordinated cellular base-station trans- 
missions |2), cooperative synchronization of pulse-coupled 
clocks via spatial averaging [16], consensus synchronization 
of clock drifts from carrier frequency measurements [17], and 
carrier-phase based pairwise ranging and synchronization [ 18 1. 
Experimental results in |[T8l show carrier-phase techniques 
can achieve offset compensation to accuracies on the order 
of 10 ps. 

This paper describes a non-hierarchical approach to ad 
hoc network synchronization based on random bidirectional 
pairwise message exchanges between nodes. Our approach is 
described in two steps: (i) drift compensation and (ii) offset 
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compensation. For clarity we present these steps as sequential 
tasks, but they can also be performed simultaneously using 
measurements at the MAC layer and/or physical layer. Since 
our approach is non-hierarchical and based on random bidi- 
rectional pairwise message exchanges, it does not require any 
coordination among the nodes and can be embedded into 
existing network traffic. Our approach is especially amenable 
to physical layer synchronization techniques which glean drift 
and offset measurements from existing data packets. All of the 
nodes in the network can run exactly the same algorithm, ir- 
respective of the network topology. Since our non-hierarchical 
synchronization algorithm is not based on local averaging, the 
overhead of our approach is potentially less than the consensus 
techniques described in Q), lfP3ll - lfT5ll . We provide analysis 
proving non-hierarchical synchronization exhibits monotonic 
expected convergence under broad conditions on the net- 
work topology and clock update stepsize. Numerical examples 
showing convergence and divergence of the proposed non- 
hierarchical network synchronization technique under different 
network topologies are also provided. 

The rest of the paper is organized as follows. We first 
introduce a probabilistic messaging model and the relevant 
local clock parameters in Section [TT] Then, we describe a 
non-hierarchical network synchronization technique based on 
random pairwise exchanges in Section [HI] Convergence results 
are also provided in this section. Numerical results are given 
in Section |IV] followed by conclusions in Section [V] Proofs 
of the theorems are given in the Appendices. 

Notation : Vectors and matrices are denoted by boldface 
letters. denotes the N x N identity matrix. The vectors 
In, Ojv, and ejy denote a vector of all ones, a vector of all 
zeros, and a vector of all zeros with a one in the last position 
respectively, all in M. N , ||-|| represents the Euclidean norm of 
the enclosed vector. We use E {•} and ( ) T for expectation and 
transposition. 

II. System Model 

We assume a time-division multiplexed network of N 
nodes with transmit/receive topology specified by a probability 
matrix P with i, j th entry pij corresponding to the probability 
that node i initiates an exchange of messages with node j 
at any particular instant in time. The case pij = pj.i = 
corresponds to the situation where node i and node j do not 
communicate. Note that p i;i = for all i = 1,...,N and 
J2iJ2jPi ,j — 1- We do not necessarily assume pij = Pj^; 
for example, pij > and pj : i = corresponds to the case 
where i initiates message exchanges with node j but node j 
never initiates message exchanges with node i. 

Since all of the channels in the system are time-division 
duplexed (TDD), we assume reciprocal propagation delays 
= ipj,i m ea ch link. Basic electromagnetic principles 
have long established that channel reciprocity holds at the 
antennas when the channel is accessed at the same frequency 
in both directions [19|. Channel reciprocity can also be quite 
accurate at intermediate-frequency (IF) and/or baseband if a 
reciprocal transceiver architecture is used [20] and can be 
further improved through transceiver calibration techniques to 
remove I/Q imbalance effects lED . Il22l . 



The nodes in the network do not possess a common notion 
of time. The following section presents a model of local and 
reference time that will be subsequently used in the description 
and analysis of the pairwise synchronization protocol. 

A. Reference Time and Local Time 

The focus of this paper is the description and analysis of 
a synchronization technique for devices in a wireless ad-hoc 
network. To support this focus, it is necessary to explicitly 
present a model of local time at each node and describe 
how the local time at each node relates to a notion of 
"reference" time. The notation t refers to the reference time, 
i.e. the "true" time, in the system. All time-based quantities 
such as propagation delays and/or frequencies are specified in 
reference time unless otherwise noted. 

None of the nodes have knowledge of the reference time t. 
The local time at node i modeled as 

U = t + Ai(t) 

where A j (t) is a non-stationary random process that captures 
the effect of clock drift, fixed local time offset, local oscillator 
phase noise, and frequency instability |23l . Over short time 
periods, a reasonable first-order model of local time can be 
written as 

U = fat + Ai 

where fa represents the nominal relative rate of the clock at 
node i with respect to the reference time and Aj is the local 
clock offset at t = 0. None of the nodes in the network have 
knowledge of fa or Aj. 

III. Non-Hierarchical Network Synchronization 

This section describes a non-hierarchical technique for 
network synchronization that allows each node in the network 
to arrive at a common clock drift fa and clock offset Aj. The 
goal is not to force fa = 1 and A,; = 0. Rather, as is often the 
case in ad hoc network synchronization [1], the goal is to drive 
the clock drifts and offsets to common values (3 and A across 
the network. For conceptual simplicity, we describe our non- 
hierarchical approach to network synchronization as a two-step 
process: (i) drift compensation and (ii) offset compensation. In 
practice, both drift and offset compensation can be performed 
simultaneously, since pairwise drift estimates can be inferred 
"for free" from normal network traffic. 

A. Step 1: Drift Compensation 

Since the reference time t is unobservable, this section 
develops a drift compensation framework in the context of 
pairwise drift estimates and compensation. In timeslot k, the 
pairwise clock drift between node i and node j is defined as 

There are N(N — 1) such pairwise drifts in the network. Since 
j3j t i[k] = —fa,j[k], we can define the network pairwise drift 
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vector at time k as 



where q t is the i column of Q N for i = 1, . . . , N with 



P[k] = 



d JV(JV-1)/2 



(1) 



where P^k] = [0i,j[ft], . . . ,/3 j -i J [k]] T e W^ 1 is a vector 
of pairwise drifts /3ij with respect to node j for all i < j. 
Furthermore, since /3jj = /J^jv — /^jv for all i,j pairs, we 
can write 

P[k]=Q N P N [k] (2) 

oN(N-l)/2xN- 



where Q N e 
A^ = 4 node case we have 



/3[fc] 
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In general, Q N has elements equal to — 1, 0, or 1. Appendix lAl 
provides a recursive definition and lists several relevant prop- 
erties of Q N . 

Since nodes derive their symbol rate and carrier frequency 
from the same local oscillator that drives the local clock, any 
message between a pair of nodes in the network allows for 
the estimation of pairwise clock drift at the physical layer 
through carrier frequency or symbol rate offset estimation. 
Pairwise clock drift can also be estimated at the MAC layer 
through observing multiple timestamped messages from an- 
other node in the network. Without restricting ourselves to 
a particular method for pairwise drift estimation, we simply 
assume that when the event occurs that nodes i and j form a 
sender/receiver pair, node i adjusts its local clock rate based 
on its local estimate of its pairwise drift with node j. More 
precisely, node i forms the estimate 0^ [k] and adjusts its local 
clock drift 

Pi[k + l]=0 i [k]+rf j , i [k] (3) 

where fi > is a stepsize parameter. All other local clock 
drifts in the network (including node j, the receiver node) are 
not updated, i.e. 

t [k + l]=0 t [k] W^i 



Note that the local drift adjustment at node i affects the 
pairwise drifts fijj and fiij for all j. Hence, we can represent 
the update to the pairwise drift vector as 



P[k+l}=p[k]+ f ib i , j [k] 



(4) 



where bij[k] denotes the pairwise drift adjustment vector 
caused by (01. As a specific example, suppose AT = 4 and 
= (2, 3). The pairwise drift adjustment vector that occurs 
when node 2 updates its local drift can be written as 

b 2 ,3 [k] = 03,2 [k] [-1, 0, 1, 0, 1, 0] T = /3 3)2 [k]q 2 (5) 



Q N = [Qn Qn] e 



sN(N-l)/2xN 



where q N = [0, . . . , 0, — 1, . . . , — 1] T is a vector of zeros with 
the last N — 1 elements equal to — 1 . 



B. Step 2: Offset Compensation 

While compensating for oscillator drift syntonizes the nodes 
in the network, it is not sufficient for synchronization because 
the fixed clock offsets among the nodes in the network are 
not corrected. To correct these offsets in a pairwise manner, 
we assume that the pairwise drift between nodes i and j is 
negligible and that node i and node j form a sender/receiver 
pair and exchange messages such that node i can disambiguate 
its pairwise clock offset with node j from the propagation 
delay i/j^j = ipjj. As one example of how this can be achieved 
0, consider the exchange of timestamped packets as shown 
in Figure Q] below. Given a packet transmitted by node i in 
local time ^ a) , it arrives at node j in local time tf = t} 



timestamps 



(d) 



The response from node j contains the local 
and tj and arrives at node i at local time 
Aj — Aj. After receiving the response, 



node i can compute the pairwise clock offset to node j as 



= A 3 - Ai = A Jti . 



After an exchange of messages with node j, node i forms the 
estimate Aj,i[k] and adjusts its local clock offset 



Ai[A + l] = Ai[k] + nA j . 



(6) 



where fi > is a stepsize parameter. All other local clock 
offsets in the network are not updated. We can represent the 
update to the pairwise offset vector A[fc] g K^f^^ 1 )/ 2 as 

A[k + 1] = A[k] + fidi tj [k] 

where dij[k) denotes the pairwise offset adjustment vector 
caused by ©. 



node i 
timebase 



node j 
timebase 



Fig. 1. 



h = t + Ai t 2 = t + A 2 

Sender/receiver two-way message exchange. 



In this context, offset compensation is conceptually identical 
to drift compensation as discussed in the previous section. The 
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pairwise offsets and corresponding update vectors have the 
same form as the pairwise drift updates in 



C. Convergence Analysis 

We use 1 1/3 [fc] HI and ||A[fc]||2 as a measure of the overall 
network pairwise drift and offset alignment at time k, re- 
spectively, with smaller values corresponding to better overall 
pairwise alignment. Since pairwise offset updates occur in 
the same manner as pairwise drift updates, our analysis here 
will focus on convergence of the pairwise drifts with the 
understanding that the theorems also apply to the convergence 
of pairwise offsets by substituting A[fc] for j3[k] and dij[k] 
for bij[k}. 

In general, pairwise drift compensation does not neces- 
sarily improve the overall pairwise drift alignment metric. 
For example, suppose f3[k] = [—1,-2,-1,1,2,3] with 
||/3[fc]||| = 20. Nodes 2 and 3 form a sender/receiver pair and 
node 2 adjusts its local drift by 183,2 [&] = —02,3[k], causing 
the drift adjustment vector b 2 ^[k] = [-1,0, 1, 0, 1, 0] T . With 
/1 = 1, the resulting pairwise drifts at time fc + 1 can be written 
as j3[k + 1] = [-2,-2,0,l,3,3] T with \\j3[k + 1]||| = 27. 
In this example, even though node 2 has perfectly corrected 
its drift with respect to node 3, the overall network drift 
alignment metric has become worse. Under certain conditions 
on p and P, however, we will show that the drift update 
prescribed in with perfect estimates $j,i[k] = /3j,i[h] = 
—Pi j [k] results in monotonic expected convergence such that 
E{||/3[fc + l]||2|/3[fc]} < ||/3[fc]||| for all 0[k] ± 0. 

Prior to stating the convergence theorems, we will define 
some relevant quantities. First, define the mapping 

3-1 

m(i,j) = i + ^{k - 1) 

fc=2 

for integer i,j satisfying 1 < i < j < N, Note that m(i,j) G 
{1, . . . , N(N — l)/2} corresponds to the position of the entry 
in /3[k] containing the pairwise drift (3ij[k]. 

In the case with perfect estimates such that $j,i[k] — 
— — Pi.j[k], the pairwise drift adjustment vector bij[k] 
is linear in /3[fc]. Hence, we can write 



where Eij G 



bij [k] = - qifcj [k] = —Q N E it j(3[k] (7) 
NxN{N-i)/2 j s a ma i- r i x f a jj zeros except 
for a sign(j — i) in position (i,m(i,j)) if i < j or position 
(i,m(j,i)) if i > j. The matrix Eij effectively selects q i 
from Qn and Pij[k] from /3[k] when i < j or —f3 iy j[k] = 
f3j t i[k] from /3[k] when i > j. 
From (O, we can write 

N N 

E{b i , j [k]\m} = EE^w 

i=i j=i 

(N N \ 
<=i 3=1 j 

= -R/3[k]. (8) 



Also note that 

Win. 



\\9iPj,i[k]\\i = {N-l)0fM] 



where we used the fact that q { has N — 1 elements equal to 
±1 and that Pj,i[k] = —0i,j[k]. It follows that 



N N 



*=1 3=1 

= (3 T {k}S(3{k] (9) 
where S G T$>n(n-i)/2xN(n-i)/2 ifJ a diagona i matrix with 
entries given as S m ^ }m ( i}j ) = (N - l){pij + Pj,i) for 1 < 
i < j < N. 

Theorem 1 (General Monotonic Expected Convergence): 
Given f3[k] ^ 0, j3j i = f3j t i, and local drift updates specified 
by ©. Then E {\\(3[k + 1]||1 1 /3[k]} < \\P[k]\\j if and only if 



Q N [R + R -fj,S)Q 



IN 



(10) 



is positive definite. 

A proof of this theorem is provided in Appendix |B] 

Theorem Q] provides an implicit condition on the stepsize /i 
sufficient for monotonic expected convergence. It is straight- 
forward to check ( fTOb numerically for a given stepsize /i. 

In networks where each node's transmission range exceeds 
the geographic span of the network, it is reasonable to model 
the network as having equiprobable sender/receiver pairs such 
that pij = jy7]^rn f° r a U * j- The following theorem 
establishes an explicit condition on the stepsize /1 for this case. 

Theorem 2 (Equiprobable Monotonic Expected Conv.): 
Given (3[k] ^ 0, fyj — local drift updates specified by 



N(N-l) 



for all 



d3j, and sender/receiver probabilities pij — 
^j. ThenE{||/3[fc+l]||||/3[fc]} < ||^[fc]||| if and only if 



H < 



N 



N - 1 



A proof of this theorem is provided in Appendix [C] The proof 
relies on several properties of Q N listed in Appendix lAl 



D. Discussion 

As a non-trivial example of a case where pairwise drift 
and pairwise offset compensation do not exhibit monotonic 
expected convergence for any /1 > 0, consider an N = 3 node 
network with 

0.9 
0.05 
0.05 

Note that all nodes initiate message exchanges in this network, 
hence each node adapts its local drift/offset at time k with 
nonzero probability. Straightforward calculations give 
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R+R , - t iS)Q 
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3.7- 1.9/x 
-0.9 



-0.9 
0.2 - O.l/i 



It can be numerically verified that one eigenvalue of this matrix 
is negative when fi = and that both eigenvalues decrease as 
fi, increases. Hence, there does not exist a stepsize p, > 
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such that this network satisfies the conditions for monotonic 
expected convergence. 

In networks that admit monotonic expected convergence for 
some stepsize /i > 0, it is possible to derive from (TT4T > in 
Appendix [B] the value of /i that provides the largest expected 
convergence step. Rewriting ( TBl i as 

g{ji) = -2fi(3 T N [*] QlRQ N (3 N [k] + ^ 2 f3 T N [k]Qj ! SQ N f3 N [k] 

the value of /.i that minimizes g(fi) is 

_ (3 T N [k]Q T N RQ N (3 N [k] 
M ° pt f3 T N [k]Q T N SQ N f3 N [k] 

which reduces to /i opt = 2 (jv^i) m equiprobable case 

since Q j^RQ jy = j^z\Q nQ n Q n^Q n = jqQ nQ n 
as shown in Appendix [C] In the non-equiprobable case, the 
optimal stepsize adapts with the current state of the pairwise 
drift vector. Exact calculation of /i opt requires knowledge of 
all of the pairwise offsets with respect to node N, hence it 
is of limited practical utility in the non-equiprobable case. 
Nevertheless, this result suggests that nodes could have a 
heuristic for adapting the local stepsize to start at a large value 
and become smaller as the network becomes more closely 
synchronized. 

Finally, although pairwise drift and offset compensation 
allow the nodes in an ad hoc network to achieve consensus 
on the drifts and offsets such that /3j — >• ft and Aj — >• A for 
all i, it is worth mentioning that the technique can also be 
used to synchronize to an external source of reference time 
if one or more nodes in the network have access to reference 
time (via, e.g. GPS). The nodes that have access to an external 
source of reference time simply do not behave as "senders" 
in the network. This forces the other nodes in the network 
to adapt to the reference time. Hence, if the conditions of 
Theorem Q] are satisfied, the network will exhibit monotonic 
expected convergence toward synchronization with an external 
reference time. 



IV. Numerical Results 

The numerical results in this section assume a network 
with N = 10 nodes, i.i.d. Gaussian distributed initial clock 
offsets Aj[0] with standard deviation 5 ms, and i.i.d Gaus- 
sian distributed initial drifts /3j[0] with standard deviation 
100 /is/iteration, for i = 1, . . . , N. In iterations k = 0, . . . , 99, 
no synchronization updates occur. During this time, the pair- 
wise drifts remain constant and the pairwise offsets tend to 
grow. For iterations k — 100,..., 499, the drift compensa- 
tion algorithm runs with randomly selected sender/receiver 
pairs with probabilities specified by P. For iterations k = 
500,501,... the offset compensation algorithm runs, also 
with randomly selected sender/receiver pairs with probabilities 
specified by P. In the following, we consider two network 
topologies: (i) a fully-connected network with equiproba- 
ble sender/receiver pairs and (ii) a partitioned network with 
equiprobable sender/receiver pairs in each partition and a 
single gateway node connecting the partitions. 



A. Equiprobable Sender/Receiver Pairs 

In this example, pij — w^rn = M ^ or a ^ ' ^ 3- 
Based on Theorem [2] both drift compensation and offset 
compensation will exhibit monotonic expected convergence if 

Figure [2] shows a single realization of the pairwise drift 
(3[k] and pairwise offset A[k] processes for k — 0, . . . , 1000 
and fixed stepsize /i = 0.5. The effect of the uncompensated 
drifts is evident in the first 100 iterations where we see linearly 
increasing pairwise offsets. In this example, the pairwise drifts 
and the pairwise offsets converge to values close to zero 
within a couple hundred iterations of when their respective 
compensation algorithms commence. 

Figure [3] shows the pairwise synchronization metrics 
and ||A[fc]|||, averaged over 1000 Monte-Carlo re- 
alizations. These results numerically demonstrate the mono- 
tonic expected convergence property of the synchronization 
algorithm for fixed stepsizes \i = 0.1,0.5,1. As expected 
from the analysis in Section IIII-DI the case when n = 0.5 
provides the fastest convergence rate since this is the closest 
stepsize value to ^ opt = 2 (N-i) ~ 0-56. The case when 
/i = 1.2 > i2 demonstrates divergence of the synchronization 
algorithm when the stepsize is too large. 




-4 1 1 1 1 1 1 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 

iteration k 



0.04 




_0.04 1 1 1 1 1 1 1 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 

iteration k 



Fig. 2. N = 10 node synchronization example with equiprobable trans- 
mit/receive pairs and fixed stepsize /i = 0.5. 

B. Partitioned Network with Gateway Node 

In this example, the N = 10 node network is partitioned 
into sets Si = {1,2,3,4,5} and S 2 = {5,6,7,8,9,10} 
where communication between members of the same set has 
probability pij = and communication between members 
of different sets has probability pij = 0. Figure [4] illustrates 
this case. Observe that node 5 is the only member of both sets, 
hence we refer to node 5 as the "gateway node" even though 
node 5 has no special functionality. From Theorem Q] we can 
numerically determine the upper bound on the stepsize to be 

< 1.11. 

Figure [5] shows a single realization of the pairwise drift 
f3[k] and pairwise offset A[fc] processes for k = 0, . . . , 1000 
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Fig. 3. Empirically averaged pairwise drift and pairwise offset norms for Fig. 5. N = 10 node synchronization example with partitioned network 
N = 10 node synchronization with equiprobable transmit/receive pairs. with gateway node and fixed stepsize [i = 0.5. 
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Fig. 4. N = 10 partitioned network topology with node 5 serving as the 
gateway node. Each edge in the graph corresponds to a sender/receiver pair 
with probability p = in either direction. 

and fixed stepsize /i = 0.5. As in the equiprobable case, 
the pairwise drifts and offsets converge to very small values 
within a few hundred iterations of the respective compensation 
algorithms, although convergence is somewhat slower in this 
example when compared to Figure [2] 

Figure [6] shows the pairwise synchronization metrics 
||/3[fc]||2 and ||A[fc]||2, averaged over 1000 Monte-Carlo re- 
alizations. We again see divergence when /i — 1.2 > 1.11 and 
mono tonic expected convergence for /i = 0.1, 0.5, 1. The main 
difference in this example with respect to the equiprobable 
case is that the convergence tends to be slower since the two 
sets only exchange messages through node 5. In particular, in 
the fj, = 0.1 case, the convergence is so slow that the drifts 
are not insignificant before offset compensation begins. 

V. Conclusions and Extensions 

This paper presents a non-hierarchical approach to time 
synchronization in wireless sensor networks that has low 
overhead and can be implemented at the physical and/or 
MAC layers. Our approach is based on random pairwise 
bidirectional messaging and allows all nodes to use exactly the 
same distributed algorithm without requiring local averaging 
of measurements from other nodes. Analytical results show 
that the non-hierarchical approach can provide monotonic 
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Fig. 6. Empirically averaged pairwise drift and pairwise offset norms for 
N = 10 node synchronization with partitioned network with gateway node. 

expected convergence of both drifts and offsets under broad 
conditions on the network topology and local clock update 
stepsize. Numerical results are also presented verifying the 
analysis under two particular network topologies. 

Potential extensions of this work include: (i) convergence 
analysis for simultaneous drift and offset compensation, (ii) the 
development of explicit bounds on the stepsize /i for network 
topologies other than equiprobable, and (iii) analysis and 
simulation of non-hierarchical synchronization with stochastic 
clocks and/or drift and offset estimations errors. 

Appendix A 
Properties of Q n 

This appendix provides a recursive definition and several 
relevant properties of the matrix Q N implicitly defined in (0. 
Let l^v G be a column vector of all ones. For N > 2, 
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Q N has the recursive definition 



Q 



N 



Qn-i Qn-i 



f N-l 



where Q 2 = 1. Recalling q N = [0, . . . , 0, -1, . . . , -1] T e 
R /v(jv-i)/2 witn the last ]\[ _i elements of q N equal to -1, 

we note that 



QnQn 



lJV-1- 



(11) 



Hence, 

QnQn 



T Qn-\Qn-1 + -fjV-2 
_lT 

NIff-i — lAT-lljV-l 



— l/V-2 

AT- 1 



(12) 



where the last equality can be shown by induction. 

Furthermore, the product QJjQn nas N — 2 eigenvalues 
equal to N and one eigenvalue equal to one. To see this, define 



G 



s (/V-l)x(/V-l) 



-T/V-2 — ljV-2 

where V is full rank. By the matrix inversion lemma [24|, we 
can write 



v- 1 



1 



N - 1 



(iV - 1)I N _ 2 - lAr_ 2 l^_ 2 

_lT 
- L /V-2 



— l/V-2 
-1 



so that indeed VV^ 1 = In-i- 

Recall that e^-i G is a column vector of zeros with 

a one in the last position. Then V(ejv-ie^_ 1 )V _1 is the 
outer product of the last column of V and the last row of 
V , hence 

VieN^el^V- 1 = 1 J—l N _ 1 ll_ 1 . 

Consequently, we can write 

QnQn = NIn—i — l/y-iljv-i 

= NVV- 1 -(N- ^VieN^eJf^V- 1 
= V (iVJiv-i ~(N- l)e A r_ ie T_ 1 ) V- 1 

where the diagonal matrix 

NI N ^ - (N - Veir-iejf^ = diag(Y, . . . , JV, 1) (13) 
contains the eigenvalues of QJjQn- 

Appendix B 
Proof of Theorem 1 

Proof: Given the pairwise drift vector f3[k] as defined 
in ([Tji and pairwise drift update vectors as defined in 

(HJi with stepsize fi. For monotonic expected convergence, we 
will show E{||/3[fc + l]||^|/3[fc]} < ||/3[fe]||| for all (3[k\ e 
range(Q^) satisfying f3[k] ^ where the expectation is taken 
over all N(N — 1) possible sender/receiver pairs. Since 



|0[* + l]|la 



\\j3[k} + fibi 



and f3[k] — Q N (3 N [k], we have monotonic expected conver- 
gence if and only if 



2(i/3jf[k]QlE{b itj [k} | 0[k]} + M 2 E 



P[k}} <o. 



From ([8} and (Q, this condition can be be equivalently 
expressed as 

P N [k]Q T N {»S - 2R) f3[k}Q N (3 N [k] < 0. (14) 

Note that R is not necessarily symmetric. Hence, (IT4b is true 
for all f3 N [k] ^ if and only if 



A = Q T N ^ 2R -V S ) + ( m -V S Y 

= Qjf{R + R T - fiS)Q N 
is positive definite. 

Appendix C 
Proof of Theorem 2 



Q 



N 



Proof: Using the notation and definitions from Ap- 
pendix |B] with the additional assumption of equiprobable 
sender/receiver pairs such that pij = p = N ^_^ for all 

we will show that 



A = Qj ! {R + R T - fiS)Q N 

is positive definite if and only if /i < j^zj- Our approach is to 
show A = f(ti)QjfQ N under the equiprobable assumption 
and to use properties of the Q N matrix from Appendix lAl to 
conclude that A is positive definite under this assumption if 
and only if /(/i) > 0. 

Observe that entry (m(i,j),i) of Q N is equal to +1 and 
entry (m(i,j),j) of Q N is equal to —1, for all 1 < i < j < 
N, with all other entries of Q N equal to zero. Hence 



N 



Q 



N 



EE*- 



where E* 



v *=i j¥* 

£ H-iVxiV(JV-i)/a j s a ma i; r i x f a jj zeros except 



for a sign(j — i) in position (i,m(i,j)) if i < j or position 
(i,m(j,i)) if i > j. Under the equiprobable assumption, R 
can be written as 

N 

E E Ei 'i 

i=l jjti 
Qn(Qn) T 

QnQn + QnQn 



1 


N(N - 


1) 


1 




N(N — 


1) 


1 




N(N — 


1) 



Note that R is symmetric in this case, hence R + R T — 2R. 
We can further compute 

2Q N RQ N = — [QJjQnQnQn + QnQnQnQn 



N(N - 1) 
2 



N(N — 1) 
(-l*r-i)(-l^_i)) 



(ATIjv-i - liv-ill_i) 



AT(AT- 1) 
_ t QnQn 



N 2 I N ^ - JVl Ar _ 1 lT / _ 1 
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where we used ( fTTT i and (fT2b from Appendix lAl in the second 
equality and ( fT2b again in the final equality. 

In the equiprobable case we also have S — 2p(N— 1) ijv-i, 
hence 



It follows that 



Q N {R + it - fiS)Q N = 



N - 1 



which, from the fact that the eigenvalues of Q N Q N are strictly 
positive as shown in [\3[ , is positive definite if and only if 
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